Effective elasticity of a flexible fllament bound to a deformable cylindrical surface 
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We use numerical simulations to show how a fully flexible filament binding to a deformable 
cylindrical surface may acquire a macroscopic persistence length and a helical conformation. This is 
a result of the nontrivial elastic response to deformations of elastic sheets. We find that the filament's 
helical pitch is completely determined by the mechanical properties of the surface, and can be tuned 
by varying the filament binding energy. We propose simple scaling arguments to understand the 
physical mechanism behind this phenomenon and present a phase diagram indicating under what 
conditions one should expect a fully flexible chain to behave as a helical semi- flexible filament. 
Finally, we discuss the implications of our results. 



A ubiquitous geometrical state filaments arrange into 
is the helix. Apart from some synthetic polymers [1 and 
biological filaments such as ds-DNA and actin filaments 
which spontaneously develop a helical conformation due 
to their inherent chemical structure [2 , helicity can also 
appear when a filament is bound to a cylindrical sur- 
face. This phenomenon can be observed across all length 
scales: from vine wrapped around trees, to DNA on car- 
bon nanotubes [4]. 

Although in several instances it is believed that what 
leads to the helicity of a filament is either a specific 
property of the filament or the specific interactions be- 
tween the filament and the underlying surface O HI |6] , 
there is evidence that semiflexible polymers binding non- 
specifically to cylindrical surfaces can spontaneously de- 
velop helical conformations. The arrangement of cellu- 
lose microfibrils in the plant cell wall [7^ is a nice exam- 
ple of it. Recently it has been suggested that the he- 
lix is the preferred conformation of semiflexible polymers 
when generically bound to the surface of an infinitely 
long cylinder, provided the cylinder's radius is sufficiently 
large [8]. 

One particular aspect of the problem that has not been 
studied and could be of great relevance, concerns the role 
of the deformability of the underlying surface. This prop- 
erty is inherent to biological materials, and the dynami- 
cal interplay between protein filaments and the soft cell 
membrane has been shown to be crucial in several biolog- 
ical processes ^. In fact, semiflexible biopolymers such 
as microtubules and actin filaments not only provide the 
cell with a highly dynamical scaffolding that regulates 
its shape, but they also mediate important extracellular 
interactions. Cell division [9^ and cell crawling [10] are 
two dramatic examples of it. 

Here, we explicitly consider the role of the surface de- 
formability and predict that new phenomenological be- 
havior arises when a filament is bound to it. We show 
how even a fully flexible filament, weakly bound to a soft 
tubular sheet, can acquire an effective, tunable helicity 
and a large bending rigidity. The physical reasons behind 
it are quite general, are applicable to arbitrary geome- 
tries, and can be understood by analyzing the nontrivial 
mechanical response of elastic sheets to local deforma- 



tions. 

Unlike fluid interfaces, the deformation of an elastic 
sheet away from its equilibrium shape comes with a bend- 
ing and stretching energy cost. It is easy to show p!T| [T2] 
that the ratio between stretching and bending for an ar- 
bitrary deformation of amplitude /i on a surface of thick- 
ness t scales as Eg/Ei) ~ (h/t)'^. Therefore, for suffi- 
ciently thin sheets, bending is the preferred mode of de- 
formation. This has a profound effect on the way elastic 
surfaces respond to deformations as the only stretch-free 
deformation involves uniaxial bending. Skin wrinkling 
under applied stress [13, 14 and stress focusing via d- 
cone formation of crumpled paper [12] are two beautiful 
examples of this phenomenon. 

Using simple scaling arguments it is possible to es- 
timate the extent of the deformation, /p, that arises 
when an indentation of amplitude h is imposed on a nar- 
row elastic sheet of width D. Following reference [12], 
the bending and stretching energies associated with this 

deformation scale as £^5 ~ i^5 {h/ D'^^ Dip and Eg ~ 

Ks {h'^ /Ip) Dip, where K^ and Kg are the bending and 
stretching constants respectively. The balance between 

the two terms gives Ip ~ Dh^ (Kg/Kij)'^ . A more fa- 
miliar form of this expression is obtained by plugging 
Kg ~ Yt and i^5 ~ Yt^ [11] (Y is the Young modulus of 
the surface) to give Ip ^ D{h/t)^ . 

This result can be generalized to a cylindrical surface 
of radius R (with D ^ R) diS long as /i <C i^, and sat- 
urates to Ip ~ R{R/t)^ for thin cylinders and/or large 
deformations |15[ [T6] . Either way, for a fixed cylindrical 
radius R and indentation /i, the extent of the deformation 
along the axis of the cylinder is set by the ratio between 
bending and stretching constants. 

In this paper we show how a fully flexible filament that 
generically binds to a deformable cylindrical surface can 
acquire a macroscopic bending rigidity and a specific in- 
trinsic curvature set by the mechanical properties of the 
surface and the extent of the deformation. The net re- 
sult is an effective semi-flexible chain that wraps around 
the cylinder with a tunable pitch. Using a combination 
of scaling arguments and numerical simulations we show 
how the characteristic length scale Ip is directly related 



to the pitch of the hehx, and we present a phase diagram 
showing the transition from a disordered (random walk) 
to the hehcal conformation of the filament as a function 
of its binding affinity to the surface. 

We model the elastic surface via a standard triangu- 
lated mesh [it]- The mesh is composed of TV = 14960 
nodes arranged to produce an initial configuration with 
perfect hexagonal tessellation. To impose surface self- 
avoidance we place hard beads in each node of the mesh. 
Any two surface beads interact via a purely repulsive 
truncated and shifted Lennard- Jones potential 
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where r is the distance between the centers of two beads, 
(J is their diameter, and e = 100/cbT. 

We enforce the surface fixed connectivity by linking 
every bead on the surface to its first neighbors via a har- 
monic spring potential 



^stretching ^sy "^BJ 



(2) 



Here Ks is the spring constant and r is the distance be- 
tween two neighboring beads, tb = 1.23cr is the equilib- 
rium bond length, and it is sufficiently short to prevent 
overlap between any two triangles on the surface even for 
moderate values of Ks. 

The bending rigidity of the elastic surface is modeled 
by a dihedral potential between adjacent triangles on the 
mesh: 



Ubending = Kb{l + COS (/)) 



(3) 



where (j) is the dihedral angle between opposite vertices of 
any two triangles sharing an edge and i^^ is the bending 
constant. 

The polymer is constructed as a "pearl necklace" with 
Njji = 20 monomers of diameter of am = 10a. Neigh- 
boring monomers are connected by harmonic springs 
as in Eq. |2] with the equilibrium bond length tm = 
1.18<Jm and spring constant of 120/cbT/(J^. Polymer self- 
avoidance is again enforced via the repulsive truncated- 
shifted Lennard-jones potential introduced in Eq. [ijwith 
a -^ am- Note that we do not associate an explicit bend- 
ing rigidity to the polymer which behaves as a simple 
self-avoiding random walk when bound to an infinitely 
rigid cylinder. 

The generic binding between polymer and surface is 
described by a Morse potential: 
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, r < lOo- 
, r > lOcr 

(4) 
where r is the center-to-center distance between a 

monomer and a surface-bead, vmb is bead-monomer con- 
tact distance tmb = 5.5cr and I^o is the binding energy. 
The interaction cutoff is set to lOa and 7 = 1.25/ a. 



We used the LAMMPS molecular dynamics pack- 
age [18] with a Nose/Hoover thermostat in the NVT 
ensemble to study the statistical behavior of the sys- 
tem. Periodic boundary conditions are imposed to make 
the cylinder effectively infinite. No difference was found 
when using the NP^T ensemble, with Pz = {z is aligned 
along the cylinder's axis). The timestep size was set to 
dt = 0.002ro (to is the dimensionless time) and each sim- 
ulation was run for a minimum of 5- 10^ steps. The radius 
of the undeformed cylinder was set to R = 14<j in all our 
simulations. 

The overall strategy of our numerical work is to per- 
form a statistical analysis of the system for different val- 
ues of Ks and I^o, and to understand how the configu- 
rational properties of the binding polymer are related to 
the elastic properties of the templating surface. 

Figure 1 shows for a particular value of the membrane 
bending rigidity the different phases of the polymer in 
terms of the binding constant I^o, which regulates the 
extent of the surface indentation /i, and the stretching 
constant Ks. We find a gas phase, an arrested phase and 
a helical phase. The behavior of the system in the limit 
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FIG. 1: Phase diagram of a fully flexible polymer binding 
to an elastic tubular surface for fixed Kb — 150/cbT. Three 
phases are shown as a function of Do and Kg - A: gas, B: helix, 
C: arrested phase. The direction of the white arrow and the 
shading in the B phase show the helical pitch increase with 
Kg. Dark regions indicate large pitch, whereas light regions 
represent low pitch. The inset shows three snapshots of the 
chain configurations in the three phases. 

of very large and very small indentations is clear. In the 
first case, the cylinder is effectively rigid and does not 
alter the behavior of the polymer which performs a self- 
avoiding random walk over its surface. We indicate this 
phases as the gas phase. In the second case, the poly- 
mer acquires non-helical conformations that differ from 
each other once simulations are repeated (under the same 
conditions) using a different initial configuration. This is 
indicative that the polymer becomes kinetically trapped 
and we take this as a signature that the system dynam- 
ics is becoming glassy. We call this phase the arrested 
phase. The most interesting behavior arises for moder- 



ate indentations, where the interplay between bending 
and stretching energies of the surface strongly affects the 
configurations of the polymer, and results in an inter- 
esting helical phase with pitch increasing monotonically 
with the membrane stretching cost. 

Local deformations caused by each monomer in the 
helical phase pair-up coherently to generate a smooth 
surface channel following the chain profile. Scaling ar- 
guments can be used to estimate the energy cost re- 
quired to form a channel along the cylinder axis and one 
around it. The total bending energy associate with the 
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axial configuration scales as E^ ~ i^5 [h/R^^ IR^ while 
that for the transversal configuration has a bending cost 

E^ r^ Kb {h/lp) Up, where / ~ cFmNm is the contour 
length of the polymer. 

As Ip is typically larger than R, Ip ~ R^h/t, the bend- 
ing energy balance favors configurations in which the 
polymer wraps around the cylinder to produce ring-like 
configurations. However, the stretching energy becomes 
negligible when the polymer is placed along the cylinder's 
axis, and grows as E^ ^ Ks {h/lp) Up when it is placed 
across the axis. The net result is that when the polymer 
is bound to a surface that is easily stretchable, i.e. suffi- 
ciently thick, it will spontaneously wrap around its axis. 
In the limit of an unstretchable, i.e. very thin surface, 
the polymer will align with the cylinder axis. The in- 
termediate regime is dominated by helical configurations 
which represent a balance between the two tendencies. 
By holding h constant and altering the relative weight 
of bending and stretching energies we can modulate the 
pitch of the helix and establish its dependence on the 
mechanical properties of the membrane. 

The angle formed between the axis of the cylinder 
and the direction of the polymer can be dimensionally 
related to the two natural length scales of the problem: 
the axial, /p, and the transversal, R 
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This functional form has the correct limiting behavior. 
In the stretching dominated regime ^ ^ 0, and in the 
bending dominated regime -^ 7t/2. It is important to 
notice that one should be able to modulate the helicity 
of the polymer by increasing its binding energy to the 
surface (i.e. h). However, for sufficiently large values of 
h the system can become kinetically trapped, or crosses 
over to the scaling behavior Ip ^ R^^'^/t^^'^ [16], which is 
independent of h. It is therefore clear how variations of 
h have a weak effect on the pitch of the polymer. 

To test our theoretical predictions, we performed a se- 
ries of numerical simulations in which we carefully inves- 
tigated the dependence of on the membrane stretching 
rigidity, and on the indentation h. The amplitude of the 
indentation, /i, is tuned by changing the strength of the 
monomer-bead attraction (binding energy) Dq, and can 
be estimated by computing the largest vertical distance 
among the surface beads underneath a given monomer. 



Fig. 2a shows how depends on Ks for fixed bending 
rigidity, i^5, and indentation, h. The line is a fit to the 
data obtained by using the inverse of the functional form 
in Eq. |5j 

Fig. 2b shows how depends on the binding energy 
1^0 which, for fixed Ks and i^5, and within the narrow 
range of values of I^o we explored, grows linearly with 
h. We repeated the calculation for two different values 
of Ks and fit the data with the inverse of Eq. |5J In both 
cases Eq. |5] appropriately describes the helicity of the 
polymer in terms of the elastic properties of the mem- 
brane. The inset of Fig. 2a shows the representative 
snapshots of the polymer conformations for different val- 
ues of Ks at constant h and K^j. Two important points 
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FIG. 2: (a) Variation of ^ as a function of Kg at fixed 
h ~ 1.17(7 and Kb = 150/cbT. The solid line indicates the 
fit to the data using the inverse of Eq.(5). The inset shows 
the representative helices form increasing values of Ks. (b) 
Variation of ^ as a function of the binding energy Do, at 
Kb — ISO/cbT, for two different values of the stretching con- 
stant: Ks = SO/cbT/ct^ and Ks = HO/cbT/ct^ The solid line 
indicates the fit to the data using the inverse of Eq.(5). 

need to be emphasized. (1) The physical origin of the 
disordered-to-helix transition of the chain can be under- 
stood in terms of the usual balance between the entropy 
of the filament and the energy penalty associated with a 
random, non-optimal distribution of indentations on the 
surface. (2) By going through the transition the filament 
acquires a large effective bending rigidity which results in 
a persistence length several times larger than the chain 
length. 

The jump in persistence length of the polymer can be 
best observed by measuring a function that accounts for 



4 



the periodic correlation between the monomers, as de- 
scribed in [T9l: 
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Here m is the number of monomers between particle i 
and j along the chain, and g{m^ i) is given by 



g{m,i) = 



(iV^-m-l)Efrr'(5M-S7j)' 



(7) 
where Sij = cos Oij is the cosine of the angle between 

bond vectors i and j, and Sij is the average over all such 

angles in the chain. Figure 3 shows G{m) for different 




FIG. 3: G{m) as calculated from Eq. (6) at Ks = IO/cbT/ct^, 
Kb — 150/cbT, for different values of Dq. The jump in per- 
sist ance length is observed around h ^ 0.5a. We indicate 
G{m) with dashed lines for h < 0.5a and we use solid lines 
for h > 0.5a. 

values of h ai Ks = 10/cbT/<j^ and K^ = ISO/cbT/ct^, and 
clearly indicates two distinct cases. For h < 0.5a the cor- 
relation between the relative location of the monomers on 
the surface is negligable, while for h > 0.5cr, G{m) shows 
perfect helical correlation of monomers over a distance 
that is larger than /. Since G{m) does not decay, it is ob- 
vious that our polymer is too short for a reliable estimate 
of the persist ance length in the helical phase. However, 
the persistence length clearly exceeds the chain length 
over several times. 

What limits the length of the polymer in our simula- 
tions is the large number of triangles required to describe 
the cylindrical surface. In fact, to avoid the artifacts due 



to the specific tesselation of the surface, monomers need 
to be significantly larger than the surface beads. We find 
that the dm = lOcr is enough for the monomers not to 
feel the underlying structure of the membrane. Interest- 
ingly, when the monomer size becomes comparable to the 
size of the surface beads, we find that the direction of the 
chain is biased along the main axes of the mesh. This is a 
reminder that below a certain length scale, the structural 
details of the underlying surface cannot be neglected. 

In conclusion, it is important to emphasize two things. 
The first is that the onset indentation amplitude h for 
helical conformation is typically just a small fraction of 
the monomer size (barely 5% in the case described above) 
which is not an unreasonable perturbation even for sim- 
ple membrane-bound proteins. The second is that al- 
though in this paper we have focused specifically on the 
problem of flexible chains on cylindrical surfaces, the na- 
ture of this phenomenon is quite general and is intrinsi- 
cally connected to the nonlinear response to deformations 
of elastic sheets. This behavior can be generalized to ar- 
bitrary geometries — we find that filaments also acquire 
very peculiar conformations when placed on spherical or 
toroidal deformable shells [20^ — and more importantly 
can be extended to any component adhering to the sur- 
face. We can anticipate [21] that elastic surfaces can be 
used to drive self-assembly of hard colloidal particles re- 
sulting in a variety of geometric patterns not unlike the 
ones observed with the filaments. Clearly the specific de- 
tails of the long-range correlations induced by the surface 
will depend on the surface topology, and on the physical 
constraints of the macromolecules adhering to it. Nev- 
ertheless, it is the interplay between the stretching and 
bending modes of the surface that will determine the ef- 
fective interactions between the components bound to it 
and the overall geometry of the aggregates. 

Our hope is that the results presented in our paper 
will stimulate experimentalists to further study the elas- 
tic and mechanical properties of elastic sheets and, in 
particular the long range correlations arising when par- 
ticles bind to it. 
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